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Abstract 

We present a functional renormalization group investigation of an Euclidean three-dimensional 
matrix Yukawa model with U (N) symmetry, which describes N = 2 Weyl fermions that effectively 
interact via a short-range repulsive interaction. This system relates to an effective low-energy 
theory of spinless electrons on the honeycomb lattice and can be seen as a simple model for 
suspended graphene. We find a continuous phase transition characterized by large anomalous 
dimensions for the fermions and composite degrees of freedom. The critical exponents define a new 
universality class distinct from Gross-Neveu type models, typically considered in this context. 
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I. INTRODUCTION 



Since the experimental realization of graphene jl]] there has been a tremendous activity 
leading to new intriguing phenomena in condensed matter physics The characteristic 
feature of graphene is the presence of the so-called Dirac points at the corners of the first 
Brillouin zone. At these special points, a linear dispersion for the low-energy excitations oc- 
curs Isj], closely resembling that of relativistic fermions. The massless relativistic dispersion 
leads to remarkable electronic properties. Prominent examples are the observation of the 



More- 



for studying long sought-after quantum 
, 0] and Zitterbewegung 3|. These phe- 



anomalous quantum Hall effect corresponding to a pseudospin— \ Berry phase 
over, graphene may serve as a simple toy mode 
relativistic effects |2| as, e. g. Klein tunneling 
nomena can be understood in the framework of non-interacting relativistic Dirac fermions 
which are realized in monolayer graphene on a substrate. However, for suspended graphene 
js| the situation is different and the system is strongly influenced by the large unscreened 



Coulomb interactions {9]. In what way the dynamics modifies the low-energy behavior of 
the excitations in graphene is an important open question. This parallels the situation in 
strongly interacting quantum field theories, as e. g. quantum chromo dynamics (QCD) where 
the interaction at low energies leads to the spontaneous breaking of chiral symmetry 10]. 
In that sense, graphene can be seen as a laboratory for strongly interacting fermions. For 
suspended graphene the essential question is if the Coulomb interactions are strong enough 
to drive the system close to an interacting fixed point. In the vicinity of a fixed point the 
system is governed only by the low-energy modes. The details of the underlying lattice 
theory are no longer relevant, and the theory drastically simplifies. There it often occurs 



that one has additional symmetries that are not present in the lattice theory ll|. Striking 
examples being the effective relativistic dispersion and the effective chiral symmetry for the 
low-energy theory. 

In this work, we consider the situation where the low-energy theory is defined in the 
vicinity of an interacting fixed point and we inquire specifically about its critical properties. 
In fact, suspended graphene may be expected to be close to a non-trivial quantum critical 



point 12J if the coupling is sufficiently strong. While in the perturbative regime the short- 
range interactions are irrelevant for the dynamics, at strong coupling this is not necessarily 
so. In fact the local four fermion interactions can be generated dynamically and may play an 



important role even for the long-range correlated system |13l-ll5l|. The role of the short-range 
repulsive interactions has been studied in the framework of the tight-binding model on the 
honeycomb lattice 16|, where a competition between the staggered density and non-trivial 
topological phases was found, depending on the strength of the interactions. Although 
both types of order are conceivable, for the case of suspended graphene, one expects a 



semimetal-Mott insulator phase transition 17H19|. where the chiral symmetry is 



broken spontaneously by a non-zero vacuum expectation value of the chiral condensate. 
This corresponds to a type of staggered density phase [iq|, Q that alternates on the two 
sublattices of the bipartite honeycomb lattice. 

Here, we specifically address the critical properties for the chiral phase transition using 



211 ] . In particular, we neglect the 



the non-perturbative functional renormalization group 
influence of the long-range Coulomb interactions and characterize the properties of the 
short-range repulsive quantum critical point. Our approach circumvents the problems of a 
purely perturbative approach close to criticality, and for the first time, allows us to follow 
the flow of this model into the broken phase. Introducing composite degrees of freedom 
for the order parameters, we show that our model has a continuous phase transition in the 
universality class of a three-dimensional matrix Yukawa model with U(2) symmetry. 

The structure of this paper is as follows: In Sec. [TT] and II III we introduce the low-energy 
theory for the relativistic quasiparticles in graphene and discuss its symmetry properties. We 
consider the behavior of fermion bilinears under the discrete parity P, charge conjugation 
C, and time-reversal T transformations and identify their role in the underlying lattice 
theory. In Sec. II VI we introduce the low-energy theory of spinless fermions on the honeycomb 
lattice which serves as a model for the dynamics of fermions interacting via a short-range 
repulsive interaction. In Sec. [V] within the framework of the functional renormalization 
group, we derive the flow equations for the partially bosonized model. Finally, in Sec. IVI| 
we characterize the properties of the continuous chiral phase transition and determine the 
critical scaling exponents. 



II. LOW-ENERGY THEORY 



On a substrate the low-energy theory of graphene is described by the free Lagrangian 

c = ir llx d^ a , (i) 
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with the linear dispersion of Dirac quasiparticles. The flavor index takes the values a = 
l,...,Nf and characterizes the physical spin of the quasiparticles. For a single layer of 
graphene, the number of Dirac fermions is Nf = 2. In this work, we take Nf = 1, which 
corresponds to a system of spinless fermions on the honeycomb lattice whose band structure 
can also be modeled by photonic crystals 22j-|24J. In the following, we will leave the value Nf 



unspecified as long as not stated otherwise. The low-energy excitations on the honeycomb 
lattice in two space dimensions are described in terms a Lagrangian in d = 3 Euclidean 
space-time dimensions, with the index /x = 0, 1,2. In Euclidean space-time the Lagrangian 



([T]) satisfies Osterwalder-Schrader reflection positivity 25], and the spinors ip^ = itfijo □ are 
not conjugate to ip, but instead define independent degrees of freedom. Furthermore, we 
use a reducible chiral representation for the fermions where the gamma matrices satisfy the 
Dirac algebra 

{ln,lu} = 26 IU ,, \i,v = 0,1,2, (2) 

and are given explicitly by 

7o = I ° ^ I , 7fc= I I , A; = 1,2, (3) 
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where <Jk, k — 1,2, 3 denote the 2x2 Pauli matrices. Apart from these matrices the Dirac 
algebra consists of the two matrices 

73= I ° 





that anticommute with all 7^, \x = 0, 1, 2, and with each other, as well as their combination 
735 = 5 [73>75]- Note, that these matrices do not appear in the Lagrangian flTJ which gives 
rise to a certain freedom to define the discrete symmetries {2^ (see Sec. IIIII) . 
In the chiral representation we take the states with definite chirality 

7 5 ^± = ±ijj± , (5) 

to define the excitations around the two distinct Dirac points K and K' = —K at opposite 
corners of the first Brillouin zone (valleys). The remaining two- components of the chiral left- 
and right-handed components essentially characterize the excitations on the two triangular 



1 For our Euclidean conventions see e. g. 
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FIG. 1. The bipartite hexagonal lattice with the two sublattices A and B indicated by full and 
open dots and the two inequivalent Dirac points K and K' at opposite corners of the first Brillouin 



zone. 



sublattices A and B of the bipartite honeycomb lattice. Both the decomposition of the 
honeycomb lattice into the two sublattices, and the two inequivalent Dirac points in the 
first Brillouin zone are illustrated in Fig. [TJ This identification of the spinor degrees of 
freedom in the chiral representation is equivalent to the one given, e. g. in 28]. 

The Lagrangian (TjQ) has a continuous U{2) chiral symmetry which is generated by the 
matrices 1, 73, 75, and 735 where we have the following transformation properties for the 
Dirac spinors: 



Ut 




i> ->■ e ie ip , 




(6) 






ip ->■ e i73 V , 


$ ->■ ^e i73e 


(7) 






%j> ->■ e i75 V , 




(8) 


Uy i5 




4; -»■ e i735 V , 


^ -»■ ^e" i7359 . 


(9) 



The U(2) chiral symmetry leads to a global U(2Nf) symmetry which plays an important 
role for the dynamics, as it constrains the possible interactions for the low-energy theory 
close to the critical point. Specifically, for the local four fermion interactions we can define 
a complete set of operators that respect the U(2Nf) flavor symmetry [19l. |29| 



fjs 



2N 



f 



2N f 



(10) 



and is furthermore invariant under the discrete parity P, charge conjugation C, and time 



reversal T symmetries. Up to Fierz transformations 
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19. 



29] there is no other operator that 



is compatible with these symmetries. The possible interactions that respect the U(2Nf) 
symmetry include a Thirring interaction term jso| and a flavor-singlet interaction which 
is of Nambu-Jona-Lasinio type 29(. Both terms fully parametrize the local short-range 
interactions on the honeycomb lattice. 



III. SYMMETRY PROPERTIES 

A possible instability triggered by the local four fermion interactions can lead to very 
different patterns of spontaneous symmetry breaking. Such an instability essentially leads 
to a non- vanishing vacuum expectation value (tjjT^vJj), where T^j is an arbitrary element 
of the Dirac algebra. To identify the physical role of the order parameters in the low- 
energy theory it is instructive to consider the properties of generic fermion bilinears under 
discrete symmetry transformations. This enables us to map the order parameters for the 
low-energy excitations onto the corresponding quantities in the microscopic lattice model. 
Such an identification is necessary, as different fermion representations may lead to very 
different interpretations in terms of the order parameters and the fermion masses that are 
generated dynamically by the interactions. It is this mapping that allows us to understand 
the properties of the possible phases and their relevance for graphene. 

A. P, C, and T 

Here, we take the parity transformation to reverse both spatial coordinates. It is given 

by 

Pip(x)P~ 1 = Pip(x) , x — (xq, — Xi, — x 2 ) , (11) 

with the unitary matrix P acting on the spinor components. Parity transformations should 
leave the kinetic term invariant, and we see that any operator of the form 

P E {70, «7i72, ^7o73, «7o75} , (12) 

will do the trick. However, these possibilities can lead to very different transformation 
properties for the fermion bilinears x^T^'ip. E. g. , in principle we could obtain a mass 
term iipip that is parity-odd. This is in contrast to the usual situation in three-dimensional 
relativistic field theories where one doubles the degrees of freedom to define a parity-even 
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mass term for the fermions 



27 



311 ] . Here, we want to find those order parameters that 



correspond to the physical excitations on the honeycomb lattice. Therefore, we define the 
discrete symmetry operations in such a way that they follow from the identification of the 
spinor components with the excitations at the Dirac points and on the distinct sublattices 
A and B. It is clear that an inversion about the center of the unit cell of the hexagonal 
lattice should exchange both the Dirac points and the sublattices. Since, in the chiral 
representation, the states with definite chirality correspond to the excitations around the 
two inequivalent Dirac points K and K' = —K, this leaves only two possibilities for the 
operator P, namely those that exchange states with opposite chirality: 7 and i^ols- In 
principle, we could choose any one of the two. We define P = 7n which yields the same 



2n 



so that the mass term 



transformation properties for the fermion bilinears ipT^ip as in 
iipip and 2^735^ are parity-even, whereas ip^fzip and ip^ip are parity-odd. For the chiral left- 
and right-handed components we have the transformation property 

P^WP-^^fi) . (13) 

Charge conjugation is defined as 

c^cr 1 = {^cf . (14) 

with the unitary operator C being any one of the following possibilities 

C e {72, 27071,^7273, 37275} • (15) 

This follows from the requirement that under charge conjugation iipj^ip —> —iip^^ip should 
hold. The question is how to constrain this set of operators. It is clear that charge conjuga- 
tion should leave the two Dirac points invariant. Furthermore it exchanges the sublattices A 
and B and it transforms particles into antiparticles. Again, we are left with two possibilities 
for the operator C: 72 and 27275. Here, we define C = 27275 where the fermion bilinears 
iipip, "ip^ip, and iipj^ip are even under charge conjugation, and ipj^ip is odd. For the chiral 
components we have the transformation 

C^CT 1 = a 2 {i, ± ) T . (16) 

We are left with the antiunitary time reversal @ 

T^T" 1 = Tip , (17) 



2 Note, that in Euclidean space time reversal simply complex conjugates c-numbers without changing the 
spatial coordinates. 
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TABLE I. Transformation properties of fermion bilinears under P, C, and T where 7^ = 
(70,-71,-72)- 

where the unitary operator 

T E {ij 2 , 7o7i, 7273, 727s} • (18) 
Time reversal changes both the momentum and spin of the quasiparticles, where we neglect 
the part of the operator that acts on the physical spin, given by some non-diagonal matrix in 
flavor space. It should exchange the two inequivalent Dirac points K and K' . Furthermore, 
time reversal will also exchange the sublattices, as particles are mapped to antiparticles. 
From these remarks it appears that we have two possibilities: 272 and 7275- We take T = 272 
for which the bilinears ii/jijj, V'73^, and ip'jsip are even, and 2^735^ is odd under time reversal. 
The transformation (fl7|) acting on the chiral components is given by 

T^T^ 1 = ±a 2 ^ T . (19) 

From these considerations it follows that in the chiral representation the mass term iipip 
is invariant separately under P, C, and T. All other mass terms break at least one of the 
discrete symmetries. The properties of the various fermion bilinears are summarized in 
Tab.H 

B. Antiunitary symmetries 

Apart from the antiunitary time reversal symmetry 

[V,i l2 K} = 0, (20) 

where K denotes the complex conjugation operator, the Euclidean Dirac operator with a 
possible mass term T> = 7^^ + m has another antiunitary symmetry which is given by 

[V, -i loll K] = . (21) 
8 
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TABLE II. Transformation properties of fermion bilinears under the antiunitary and discrete chiral 
transformations. 

In terms of the components of the Dirac spinor this symmetry exchanges the excitations 
on the two sublattices A and B, as well as the physical spin of the quasiparticles. As for 
the physical time reversal symmetry, we neglect the part acting on the spin components. 
Eq. ( 12~TT) can be seen as a time reversal operator acting separately at the two inequivalent 



Dirac points j(3]. While the time reversal operator satisfies (i^K) 2 = 1 and therefore defines 
an orthogonal symmetry, we have (— 27071-K') 2 = —1 which is a symplectic symmetry. For 
the low-energy theory ([I]) the orthogonal time reversal symmetry follows simply from the 
fact that we have an even number of chiral fermions in our set-up. Stated differently, this 



is just the celebrated Nielsen-Ninomiya theorem [32]. One may easily verify that both the 
orthogonal and symplectic symmetries can be constructed from two antiunitary symplectic 
symmetries acting at each of the Dirac points independently. Thus, the low-energy effective 
theory of monolayer graphene at the charge neutral point can be identified with two degen- 
erate copies of the All (Wigner-Dyson symplectic) symmetry class 33[ following the Cartan 



classification of symmetric spaces 



34) 



The commutator of both antiunitary operators T and S vanishes [T, S] = and so their 
product 

TS = 2 735 , (22) 

gives a well-defined unitary operator, exchanging only the two Dirac points. Clearly, if both 
T and S are symmetries of the theory then the discrete chiral transformation TS also defines 
a symmetry operation. In Tab. [XT] we have collected the transformation properties of the 
fermion bilinears under the antiunitary and discrete chiral transformations. 

Typically, in QCD-like theories the antiunitary symmetry of the Dirac operator is related 
to the (pseudo)reality of the fermion color representation. For these theories, where the 



fermions are massless, the antiunitary symmetry is responsible for an enlargement of the 
SU(Nf) x SU(Nj) x U(l) chiral and baryon number symmetries to a global SU(2Nf) flavor 



symmetry 



35 



361 ] . However, here, we have an even number of massless chiral fermions in a 



complex representation. Both the (orthogonal) antiunitary symmetry of the Dirac operator 
and the extended flavor symmetry are a simple consequence of the reducible representation 
for the fermions in three space-time dimensions. 



C. Order parameters 

From the above discussion it follows that the expectation value i(ipip) is invariant under 
P, C, and T whereas both (^73^?) and ('tjrygijj) are parity-odd (compare Tab. [I]). All three 
order parameters break the extended U(2Nf) flavor symmetry, generating a dynamical mass 
for the fermions. The symmetry breaking pattern is given by 37] 



U(2N f ) -» U(Nf) x U(N f ) . 



(23) 



Thus, it is clear how to identify these order parameters with the excitations in the under- 
lying lattice model. They can be mapped onto the staggered density phase where one has 



19 



20 



38 



and a bond 



an alternating density on the two different sublattices A and B 
ordered phase that corresponds to a hopping texture on the nearest-neighbor links in the 
language of the tight-binding model 39(. Indeed, a non- vanishing staggered density on the 
two sublattices breaks parity and we may associate the order parameter (^73^) with such a 
phase. The corresponding bilinear measures the imbalance in the local densities of the two 
sublattices and therefore does not mix the chiral modes. However, the order parameter for 



the bond ordered phase should couple excitations 



K and K' and on the two sublattices A and B 



28 



Doth at the two inequivalent Dirac points 



391 ] , which is accomplished by the bilinear 



i( , 0(cos a + 75 sin a)ip). The parameter a controls the hopping texture, and depending on its 
value, one obtains either a parity-conserving, or a parity-breaking type of order (see Tab. [I]). 

The condensate 2(^735^) i s m a sense special, as it leaves the chiral symmetry intact, 
however, breaks the time reversal symmetry. It corresponds to a topologically non-trivial 
ohase that relates to counter-propagating currents on the two different types of sublattices 



401 ] . This identification of the order parameters is equivalent to the one proposed in 28 ] 



where also a chiral representation for the Dirac algebra was used. A complete classification 



10 



of all possible bilinears is given in 4l|, |42 ] . 



IV. A SIMPLE MODEL: SPINLESS FERMIONS ON THE HONEYCOMB LAT- 
TICE 

Here, we focus on the case of a single flavor Dirac Lagrangian with a local four fermion 
interaction. This essentially corresponds to a low-energy theory of spinless fermions on the 



honeycomb lattice with a nearest-neighbor and next-to- nearest neighbor coupling [16]. In the 
purely fermionic description, one may expect that above some critical value for the coupling, 
the short-range repulsive interactions give rise to an instability. Here, we go beyond a mean 
field approach by introducing a set of composite fields that allow us to follow the system 
into the ordered phase. Even for this simple model there is a complex phase diagram with 
different types of order: Depending on the strength of the interactions there is a competition 
between a staggered density phase and a topologically non-trivial phase [if]]. The global 
symmetry group for this model is U(2) and the symmetry breaking pattern for the chiral 
transition is given by SU(2) — > U(l) while the topological phase transition leaves the U(2) 
symmetry intact but breaks T symmetry. 

In terms of the microscopic low-energy theory (flOj) the relevant physics for the SU (2) — >■ 
U(l) chiral phase transition is fully characterized by the fluctuations in the flavor-singlet 
channel since, for a small number of fermion species, the vector-like interactions are irrelevant 



at the transition 29]. Thus, in the following, we neglect the Thirring interaction and focus 
exclusively on the flavor- singlet term. The scalar singlet term can be brought to a flavor- 
nondiagonal form by virtue of a Fierz transformation. This defines a generalized Nambu- 
Jona-Lasinio model which serves as the appropriate low-energy description of a system of 
spinless fermions on the honeycomb lattice. To see how the repulsive interactions drive 
the system into the broken phase we integrate out the flavor non-diagonal four fermion 
interaction via a Hubbard- Stratonovich transformation. In this way we obtain a matrix 
Yukawa model with a U(N) symmetry for N = 2 species of massless, two-component Weyl 
fermions ip a and ip a , a = 1, . . . , N. For the spinless fermions on the honeycomb lattice N = 2, 
which corresponds to a single Dirac fermion Nf = 1 in the reducible representation, as e. g. 



modeled in microwave photonic crystals 



24j . The Weyl fermions couple to a Hermitian 
11 



matrix field § a b and the action of this model is given by □ 

S , [$,^,V] = J |^# a + h^ a i<$> ab ^j b + ^m 2 tr$ 2 | , (24) 

where the trace tr(- • • ) acts on the indices of the matrix field. Here, and in the following, 
we define <j) = a^d^ which belongs to the irreducible representation 7 = <t 3 , and 7^ = a k) 
k = 1,2. Eq. ([21]) constitutes the starting point for our investigation in the framework of 
the functional renormalization group. 



V. FUNCTIONAL RENORMALIZATION GROUP 

The full information about the quantum dynamics of a theory is given by the quantum 
effective action T, which is the generating functional for one-particle irreducible correlation 
functions. The functional renormalization group is a non-perturbative approach to determine 
the quantum effective action, taking into account all quantum fluctuations. Implementing 
Wilson's renormalization group idea Q|, the fluctuations are included successively by inte- 
grating out the higher modes. Thereby, one obtains the effective average action Tk with all 
the fluctuations included above the characteristic momentum scale k. The scale-dependence 
is implemented by an infrared regulator that gives a mass to the low-momentum modes in 
the theory. In the limit k — > 0, when all quantum fluctuations are included, the functional 
renormalization group yields the full effective action V. 



The flow equation for the effective average action is given by |2l|, |43] 

d t r k [ x ] = \stt \d t n k (rf 1 ^] + 7^)^} , (25) 

where t = ln/c/A defines the scale parameter, and A is some appropriate ultraviolet scale 
where we impose the microscopic dynamics. The supertrace STr in fl25l) denotes a summation 
over fields and possible internal indices, as well as an integration over momentum, while 
providing a minus sign in the fermionic sector. The second functional derivatives of the 
effective average action I 1 ! 1 define the fluctuation matrix. In the momentum representation, 



A' 

we have 



7 „ r , V 



[x] ) {p-i q) = t~yi \ r/c[x] 7 ' (26) 



We use the notation J x = J d d x for the space and j = ? d for the momentum integrals. 
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where the indices /, J label the different components of the auxiliary field x that contains 
the complete field content of our model: the matrix field $ a & and N species of Weyl fermions 
ifj a and ip a , i. e. Q 

x T (-p) = ($ T (-p) , i> T (-p) , <Kp)) , (27) 

and 



(28) 



where we have suppressed the flavor indices. Together with the infrared regulator TZ k it 
represents the full regularized inverse propagator fr^. 1 ' 1 ^ + TZkj at the scale k. 
The regulator function TZ k takes the following form 



RbM o o 
o o r f 4p) 

\ Rl k (p) ) 
where in the bosonic and fermionic sector, we have 

RbAp) = Z B,kP 2 r B ,k(p 2 ) , 
RF,k(p) = Z F)k f>r F ^{p 2 ) . 



(29) 



(30) 
(31) 



Both are described by the regulator shape functions r B)k and that characterize the 
scheme- dependence of the renormalization procedure. Since they depend only on the di- 
mensionless ratio y = p 2 /k 2 we will drop the index k in the following. The fermionic 
regulator shape function r F is taken to satisfy the constraint p 2 (l + r^) = p 2 (l + r F ) 2 and 
thus, is completely determined by the choice of tb- 

In practice, to solve (121)1) one is bound to rely on approximations for the effective average 
action Y k where one truncates the set of possible operators following an expansion scheme, 
e. g. in powers of derivatives j^] . However, such an approximation also induces a spurious 
dependence on the regulator for the full quantum effective action when the scale k is sent to 
zero 45|, |46| . In this work we therefore employ two different types of infrared regulators to 



test the regulator scheme-dependence of our results in the physical limit. We consider the 



Since the matrix- valued field $ is Hermitian, $ and $* arc not independent degrees of freedom. 
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optimized regulator [45] 



Mv)= Q-lW-y) , (32) 
and also an exponential-type regulator 

r B (y) = (exp(y) - l)' 1 . (33) 
A. Effective average action 

Our ansatz for the effective average action is given by 

T fc [$, V] = / {^,* + ^fl,fc tr (<9 M $) 2 + ^$ ab ^ + } • (34) 

In contrast to the microscopic model, we have a kinetic term for the composite field and a 
scale- dependent wavefunction renormalization both for the fermions and the bosons. Thus, 
we include the bosonic fluctuations that give a non-trivial momentum structure for the 



fermion interactions jiol . 



~^\ . Clearly, our ansatz (134"]) goes far beyond a simple mean-field 



approximation where one neglects the contributions from the composite degrees of freedom. 

For the U(N) matrix-model the effective average potential £4($) is a function of the 
invariants of the U(N) symmetry group. For system of N = 2 Weyl fermions there are 
exactly two invariants, one that is linear = tr$, and another one that is quadratic 
pk = | tr $ 2 in the fields. A non- vanishing vacuum expectation value for the composite 
field $ a f, signals the dynamical generation of a mass for the fermions. Depending on the 
flavor structure of the matrix field we have different types of order: The chirally broken 
phase corresponds to a vacuum configuration which is either non-diagonal, or diagonal non- 
uniform (so that the trace vanishes, i.e. tr$ = 0). On the other hand, $ j ~ 5 a b in the 
non-trivial topological phase which we will not consider here. Thus, to investigate the nature 
of the chiral phase transition we may neglect the fluctuations of the cx^-field. However, we 
want to emphasize that this approximation is no way essential for the following calculations. 

We expand the effective average potential in powers of pk around the minimum p ,fc> given 
by the scale-dependent vacuum expectation value: 



' m a ,r 



U k (pk) = ml(pk - po,k) + ^2 ~ ■ ( 35 ) 



n=2 



This approximation captures all the relevant fluctuation at the chiral phase transition. In 
the symmetric regime the vacuum expectation value po,fc is zero whereas in the chirally 

14 



broken phase po,/t 7^ and the mass becomes zero. In our ansatz we include the first set 
of irrelevant operators according to a naive power counting with respect to the canonical 
mass dimension. 

B. Flow equation for the effective average potential 

To extract the flow equations for the parameters and couplings in the effective average 
action (|34p one has to project the functional flow given by the r. h. s. of 025|) onto the 
corresponding operators. For the couplings that appear in the effective average potential 
this is done by evaluating the second functional derivative rjj. 1 ' 1 ^ in a constant background 
configuration of the matrix field $ a &. 

The flow equation for the effective average potential receives contributions both from the 
boson and fermion degrees of freedom 



To keep the notation clear, we will drop the fc-index in the following. Where necessary, we 
will revert to our original notation. 

From (j25l) we obtain the boson contribution to the effective average potential 



where the full regularized boson propagator Pb is given in the Appendix |A] We sum over all 
mass eigenvalues M Bi of the mass matrix as given by the second derivatives of the potential 



In general this expression will be momentum- dependent. However, for the calculation of 
the effective average potential with a constant background configuration $ a b all momentum- 
dependence in (1381) drops out. The mass eigenvalues are given in Appendix ICl 

For illustrational purposes we perform the derivation of the flow equation in the fermion 
sector explicitly by integrating out the fermions in the action. To evaluate the fermion 
contribution to the effective average potential it is useful to write the action in terms of the 
four- component spinors 



d t U k ($)=d t U B>k ($)+d t U Fik ($) . 



(36) 




(37) 




(38) 




(39) 
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which is constructed from the two independent degrees of freedom ip and ip. In that case, 
the fermion bilinear part of the action takes the form 

f m T -(- q )V ab {q)m\q) , (40) 

in momentum space. Note, that the domain of integration is restricted to positive frequencies 
go > to counteract the doubling of degrees of freedom that comes from switching to the 
four-component spinors ( I39p . In this basis, the inverse regularized fermion propagator is 
given by 

V ab (q) = Z F (l + r F )[ & )s ai + ih(q)[ , (41) 

where h(q) = h(—q, q) denotes the momentum-dependent Yukawa couplingjfl As for the 
boson contribution, we evaluate the inverse propagator T> ab (q) in a constant background 
field Q ab . Performing the integration over the Grassmann fields, the fermion contribution to 
the potential takes the form 



U F = - ]ndetV(q) , (42) 

Jq,(qo>0) 

where the determinant acts on the flavor and spinor indices. To evaluate this expression we 
put T> a b(q) in standard diagonal form. That is, by a unitary transformation we diagonalize 
$ a 6 = $<A&, so that 

V ab (q) = V a (q)5 ab , (43) 
and the determinant in ( 142]) can be written as 

detP a = (Z 2 F {\ + r F fq 2 + %) 2 <^) 2 . (44) 

With this result, the fermion contribution to the flow equation ( )36|) becomes 

d t U F = -2 j q 2 Z F (l + r F )d t (Z F r F ) ^P F (M Fa ) , (45) 

and one may easily verify that this is just the result that is obtained when the supertrace 
in (I25p is computed directly, using the definition of the full regularized propagators given in 
Appendix [A] and the regulator R F . Here, the masses M Fa denote the iV eigenvalues of the 
N x N matrix h(q) 2 $ 2 . 



5 For the momentum dependent Yukawa coupling h(—p, q) the momenta p and q denote the incoming 
fermion momenta at the Yukawa vertex. 
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Eq. ( 137)) and (145)) give the full contribution to the effective average potential. To in- 
vestigate the critical properties at the phase transition, however, it is convenient to bring 
the flow equations to a form where one may easily identify possible fixed point solutions. 
For that purpose, we switch to dimensionless renormalized quantities p = k 2 ~ d Z B p and 
u(p) = k~ d Uk(p)- Then, the flow equation for the potential is given by 



-da tt I , fj n , „ \„u-d®Uk 



d t u = -du + k~ d d t U k + (d- 2 + 7] B )pk^ 



dp 



(46) 



where we have introduced the scalar anomalous dimension r\ B = —d t lnZ B . Substituting 
our previous result this finally gives 

t u = -du + (d-2 + VB )pu' + 2v d |(iV 2 - 2)4 S V; Vb) + li B) (u' + 2pu"; r) B ) 

-2NlfU^ P h*-r lF \y (47) 



where the prime u' = ^ denotes differentiation with respect to the dimensionless renormal- 
ized field p, and Vd = (2 d+1 7r d ^ 2 T (2/<i)) . Here, we have introduced the threshold functions 
Iq and /g that parametrize the bosonic and fermion one-loop integrals contributing to the 
effective average potential. They are defined in Appendix [B] where their form is given ex- 
plicitly for the optimized regulator (1321) . The corresponding expressions for the exponential 

2l| . As stated in the previous section, the threshold 



regulator (j33j) can be found in, e. g. 



10 



functions carry the full scheme- dependence of the renormalization group equations. In that 
sense, the flow equations are universal - only the dimensionality and symmetries determine 
the flow and the regulator-dependence resides solely in the threshold functions. 

In the symmetric regime, we may derive the flow equations for the dimensionless renor- 
malized couplings e = k~ 2 Z^ l 7fi 2 and A n = k'* n ~ 1 ^ d ~ 2n Zg n \ n , n = 2, . . . ,n max , from (]4"7j) 
simply by differentiating with respect to field p, that is, we have e = u' for the mass parame- 
ter and X n = xir' for the couplings. The derivatives of the threshold functions are evaluated 
as 

^i B) (w; VB ) = - (n + 5 n ,o) Vb) , (48) 

and equivalently for ln(w; rjp)- Here, we give the flow equations for the mass parameter e, 
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and the couplings A2, and A3 in the symmetric phase: 



d t e = (-2 + VB )e - 2v d { (N 2 + l) A 2 l[ B) (e; Vb ) - 4h 2 l[ F) (0; Vf ) } , (49) 
d t X 2 = {d - 4 + 2 VB )\ 2 + 2v d { (iV 2 + 7) A 2 4 B) (e; r/ B ) - (iV 2 + 3) A 3 l[ B) (e; r/ B ) 

- ^ 4 4 F) (0;^)} , (50) 

ft A 3 = (2d - 6 + 3r/ B )A 3 - 2v d {2 (iV 2 + 25) X 3 2 1 { 3 B) (e; r? B ) - 3 (N 2 + 13) A 2 A 3 1 {B) (e; t? b ) 

+ (^ + 5)^^(6;^)-^^ (0;^)} . (51) 

The flow equations for the higher order couplings can be obtained by a simple differentiation 
with respect to the field, and are not given here explicitly 

As the system goes over into the broken phase, the scale-dependent mass parameter e 
goes to zero, and the field assumes a non-vanishing expectation value po 7^ 0, defined by 
u'(po) = 0. Due to the scale-dependence of po, which is given by 

dtpo = —r-dtu'(po) , (52) 
^2 

we get an additional contribution to the flow (147]) in the broken phase: 

d t X n = d t X n \ P0 + Xn+idtPo ■ (53) 
There, the flow equations for p , and the couplings A 2 , and A 3 are given by 

d tPo = (2 - d - V b)po + 2v d | (N 2 - 2) l[ B) (0; rj B ) + (3 + ^) l[ B) (2p A 2 ; V b) 

-^ 2 ^(^ 2; ^)}' (54) 
d t X 2 = (d - 4 + 2 VB )X 2 + 2^1 (iV 2 - 2) A 2 4 B) (0; r) B ) 

+ (3A 2 + 2 Po A 3 ) 2 l [ 2 B) (2p A 2 ; Vb) ~ ( 2A 3 + 2 Po A 4 - ^) lf ] (2p A 2 ; V b) 



X 2 

- ^ 4 F) ( §M 2 ; ^) - 4^ 2 4 F) (|m 2 ; ^ 1 , (55) 
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d t X 3 = (2d - 6 + 3 VB )\ 3 - 2vJ (N 2 - 2) (2\ 3 2 /f } (0; Vb ) - 3A 2 A 3 (0; r, B ) 



+ 2(3A 2 + 2p A3) 3 4 B) (2p A 2 ;r/ B ) 

- 3 (3A 2 + 2p A 3 ) (5A 3 + 2 Po A 4 ) l[ B) (2p A 2 ; V b) 



+ 4A 4 + 2p A £ 



2poA 3 A^ 
A, 



l[ B) (2p \ 2 ; VB ) 
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A 



A, 



(56) 



The flow equations for the higher order couplings can easily be obtained via ( 1531) . 

Recall that the relevant symmetry breaking pattern for the simple model (Sec. HVj) is 
given by SU(2) — > U(l). This is in direct correspondence to the 0(3) — > 0(2) transition 



in the three-component vector model 



n 



48l . l49j . Building on our previous remark concerning 
the universality of the renormalization group flow, we observe that in the bosonic sector 
(neglecting the fermion contributions) the flow equations for the effective potential are iden- 
tical to the flow equations for the three-dimensional 0(3) vector model 49] . We want to 
emphasize that this is a simplification that occurs only for the special case where N = 2. 
In the general case, the flow equations for the effective potential correspond to the matrix 
Yukawa model with U(N) symmetry. 



C. Boson anomalous dimension 

For the computation of the boson anomalous dimension r] B = —dt In Z B we first evaluate 
the flow equation fl25|) in a spatially varying field configuration $(a;). This is necessary for 
the projection onto the kinetic term and the wavefunction renormalization Z B . In partic- 
ular, we consider a distortion around the non-diagonal vacuum configuration = $o^a& 
characterized by a non- vanishing momentum Q, i.e. 

$ ab (x) = $ S ab + (5§e- iQx + c.c.) A ab , (57) 

V v ' 

= A(aO 

where the Hermitian matrices E and A satisfy the properties S T = — E and A T = A. Clearly, 
the fluctuations in the A-direction are orthogonal to the ground state orientation. Though 
we take only one of the possible orthogonal directions for the fluctuations into account this 
still yields a complete description of those contributions coming from the Goldstone modes. 

19 



Of course, in the broken phase, we also have fluctuations ~ A'(x)E a 6 from the radial mode 
that give additional contributions to the boson anomalous dimension. 
In momentum space the configuration (157)) reads 

$ ab (p) = $ 5(p, 0)S Qb + A(p, Q)K ab , (58) 

where we define 5(p,q) = (2ir) d 5^(p — q), and the amplitude is given by A(p, Q) = 
(6&5(p, Q) + 5&*5(—p, Q)). Taking the ansatz (1571) one may easily verify that 



d t ZB = — lim 



lim ^„ , d t T k 



(59) 



NQ^odQ 2 [a*->o 

gives us the flow equation for the momentum-independent part of the wavefunction renor- 
malization Zg. Here, we neglect all momentum-dependence of the wavefunction renormaliza- 
tion. Eq. ( 159)) defines a projection of the flow equation (1251) onto the flow of the wavefunction 
renormalization, i.e. dtT k \ Zg = dtZs- In the following, we will use this notation frequently. 

To evaluate (159)) we make use of a series expansion of the flow equation. Using the 
decomposition of the full inverse regularized propagator 

if • *> +7Z k = + JF h , (60) 

which is written in terms of T\~ containing all field-dependent fluctuations around the back- 
ground field configuration, and the background- field propagator V^ 1 that carries the explicit 
regulator dependence, we may write the flow equation (125)1 as a series expansion in powers 
of the fields: 

d t T k = l - STr d t K k V k + l - STr d t (V k T k ) - \ STr d t {V k T k f + 0{ H) ■ (61) 

Here we have defined the formal derivative operator 

d 

d t = dtKk—, yr- . (62) 

In terms of (157)) the leading contribution to the fluctuation is J 7 ^ ~ 0(A n ) and for the 
calculation of the anomalous dimension only the second order term S^T k in the fluctuation 
T k is important. Clearly, the lowest order term will not contribute, as it is independent of 
5$ and <5$* and thus yields a vanishing contribution to ( 1591 . The 5^T k term does include 
combinations of the type oc 5$<5$*, however, they are independent of momentum Q 2 . Thus, 
the lowest order term relevant in (1591 is the second order term 5^T k in the series expansion, 
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and specifically, we will need those terms in that are of linear order in the amplitude 
A. Apart from the bosonic background $ a f, we need to specify the fermionic background 
configuration, where we take 

i> = i> = . (63) 

Then, the matrix of second functional derivatives rjj. 1 ' 1 ^ becomes block-diagonal in the bo- 
son and fermion subspaces and we may treat the boson and fermion contributions to (159j) 
separately. 

We start with the bosonic sector. For matrix-valued fields, it is convenient to work in 



the non-diagonal basis for the propagators and fluctuations 50j. We follow the outline 



given above, and evaluate the second order term S^Tb = —jSTrd t (PbFb) in the series 
expansion where the index B denotes the corresponding quantities in the bosonic sector 
where we have dropped the k-mdex for clarity, i.e. Ys,k = T_b, Pb,u = Pb, etc. For that 
we need the boson propagator in the constant background configuration $ a6 = §oTj ab which 
takes the following form 

(PB) ab , cd (p) = ^ (s ac 6 bd - A{p) B +NB ^ d ) , (64) 
where we have introduced the quantities 

A(p) = Z B (l + r B )p 2 + ^, B = P ^. (65) 



Furthermore, we need the contribution from the fluctuations to linear order in A, which is 
given by 

(i^Ucdfo q) = A(p - q, Q) (S afe Af d + A ab ^ d ) + 0(A 2 ) . (66) 

With these results we can immediately compute the trace in 5^'Tb and evaluate the pro- 
jection onto the kinetic term: 

*»r B L = -n$1 ( m'to* f 4 ' , A J ■ (67) 



1zb ~ °\dp 2 J Q-+*dQ 2 J p t \A(p)A(p + Q)+NB 
This expression is proportional to the vacuum amplitude $o an d thus, the boson contribution 
to the wavefunction renormalization vanishes in the symmetric regime. 

We evaluate the contribution to the flow equivalently in the fermion subspace. For the 
non-diagonal background configuration (IBTI) the fermion propagator is given by 

/ (P i+) ) 
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where the components 

{P F ±] )Jp) = Pf(p) (Ml + r F )i>8 ab =F ih(p)d Kb) , (69) 

and 

P F ip) = (Z\{\ + r F )V + Hpf^ty 1 . (70) 
In the fermion subspace the fluctuations take the form 



-A T ab 



(F F ) ab (p, q) = ih(-p, q)A(p -q, Q)\ a ° 5(p, q) + 0(A 2 ) . (71) 

Aab 

Going through the same steps as above, that is, computing the trace in S^T F , and evaluating 
the projection (loTJj) we obtain 

6^T f \ Zb =-2 hm^ J [h(-p,p + Q)] 2 d t {z F (p)(l + r F (p))P F (p) 
x Z F (p + Q)(1 + r F (p + Q)) P F (p + Q) 

+ P h{p)h{p + Q)P F (p)P F (p + Q)} , (72) 

where p,p + Q)l 2 = h(—p,p + Q)h(—p — Q,p). Putting both results ( |67j) and (J721) 
together c^Zb = + 5^r^| z and using the definition t\b = —dtlnZs we obtain 

the evolution equation for the boson anomalous dimension: 

T] B = 16^ ^poXlm { 2 B j(e,e + 2p X 2 ]r] B ) + h 2 m{ F) ^p /i 2 ;?7 F ^ 

+ ^M 4 mf ) (|M 2 ;^)}. (73) 

( B) ( F) ( F) 

Here we have introduced the threshold functions m\ 2 , m\ and m\ that define the one- 
loop contribution appearing in the calculation of the wavefunction renormalization. They 
are given explicitly in Appendix [B] 

In the broken phase the wavefunction renormalization will get additional contributions 

49j . In principle, these terms should be taken into account, 



from the radial mode 
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however, since most of the running takes place in the symmetric regime, we do not expect 
them to play any major role. 



D. Fermion anomalous dimension and Yukawa coupling 



The derivation of the flow equation for the fermion anomalous dimension t] F = —d t In Z F 
and the Yukawa coupling h 2 proceeds in the same way as explained in Sec. IV CI Here, the 
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only difference is, that we need to choose a non- homogeneous configuration for the fermion 
fields, where in the momentum-representation we have 

V>(?)=#(<Z,Q) , $(q) = '$6(q,Q) . (74) 

The matrix field is evaluated in the constant background configuration <& ab = $o£ a f>, where 
the boson and fermion propagators are given by (j64"|) and fl69|) . respectively. 

Starting from the ansatz for the fermions (1741) we evaluate the fluctuation matrix 

ih(p-q,q)$ b (q-p)6 acl -ih(q, -q + p) ^{-q + p)S bc > 

-ih(-q + p,-p)i>l,(q - p)S b > c 

ih(-p,p- q)4> b >(-q + p)5 a , c 



V 



where we indicate the flavor indices explicitly. Recall, that the fluctuation is defined as the 
field-dependent part of the second functional derivative of the effective average action. The 
functional derivative from the left is taken with respect to the fields V'J > V'c) which de- 
fines the row indices. The primed column indices are defined equivalently via the right-hand 
functional derivative. Together with the background-field propagator Vk (see Appendix [Aj 
we evaluate the second order contribution 5T^ = — | STr d t {VkFk) 2 in the series expansion. 
A short calculation yields 

5r k ] = lj d t {[H P - q ,q)]\p B ) ab , c Ap)Mq-p){P i / ) )Jq)Mq-p) 
+ [Hq,p-q)]\PB) ab4 Ap)^(-q + p)(P { fY cb (q)^(p-q)j . 

(76) 

Inserting the expressions the boson and fermion propagators ( 16~4"|) and fIBT?]) this can be written 
in the form: 



5T 



(2) 



1 

N 



+ ih(Q)r® z ab 



[h(Q,p)Y Pf(p) 
2 



Z F (l + r F )$ a j> 
1 



N 2 - 2 



1 



A(p-Q) A(p-Q) 



NB 



(77) 



A(p-Q) A(p-Q) + NB 
Projecting this equation onto the corresponding operators in the ansatz for the effective 

average action we obtain the evolution equation for the fermion anomalous dimension 



Vd 



N d 



= Tf-jti 2 <! (N 2 - 2)m[ F 2 B) [ -^Poh' 2 , e; rj F , r] B 



N r 



AFB) 



(78) 



\ 



/ 



(75) 
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and the flow equation for the momentum-independent part of the Yukawa coupling 



d t h 2 = (d-4 + 2 VF + VB )h 2 - -h% d 1 2 ^ 

- l[ F i B) (^Poh 2 ,e + 2p X 2 ;r] F ,r] B 

- (2A 2 + 2p A 3 ) l[ F 2 B) (j^Poh 2 , e + 2p A 2 ; Vf, Vb 

32 fi I JFB) I 2 9 l 

+ ^Po« w d <^ 2 Z^ I — p ft , e; p F , p B 1 

- (y^ 2 , e + 2p A 2 ; rj F , r, B j } . (79) 

Here, we have defined the threshold functions rnf 2 B \ lv2^ ■> an d that are gi ven 

explicitly in Appendix [B] Together with the flow equations for the parameters and couplings 
of the effective average potential and the anomalous dimension for the bosons they constitute 
a closed set of differential equations that can be solved numerically. 



VI. RESULTS FOR THE QUANTUM PHASE TRANSITION 

We evolve the flow equations starting from an appropriate ultraviolet scale A to the phys- 
ical limit k —> 0. The following results have been obtained for the set of initial conditions: 
Z BA = 10~ 10 , Z FA = 1, h 2 A = A, where the mass at the ultraviolet scale m\ is taken free 
parameter, and all higher order couplings are set to zero, i. e. A n] A = 0, n = 2, . . . , n max . By 
varying the dimensionless mass at the ultraviolet scale we may tune the system across a 
quantum phase transition. Here, 8e = \e\ — eA,( cr )| measures the deviation of the parameter 
e\ from its critical value £A,(cr)- Close to 6a (cr) we find a fixed point solution for the di- 
mensionless renormalized parameters and couplings. That is, the parameters and couplings 
stay nearly constant over a wide range of scales as illustrated in Fig. |2j This is a clear 
indication for the presence of a continuous phase transition where the system displays a 
universal scaling behavior. For 5e > 0, where the mass parameter e\ is above the critical 
value 6a, (cr), the solution stays in the symmetric regime. However, starting just below e\^ cr ) 
the scale-dependent mass eventually becomes negative which signals the transition into the 
broken phase. This result does not depend on the special choice of initial conditions, that is, 
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FIG. 2. Renormalization group flow (a) for the dimensionless renormalized parameter e, the 
(rescaled) dimensionless renormalized couplings 10 _1 A2, and 10 -2 A3, and (b) for the boson and 
fermion anomalous dimensions t]b and rjp as a function of the scale parameter t = In k /A close to 
the critical trajectory. The full and dashed curves refer to initial conditions just above and below 
the critical parameter e^Ycr)) respectively. 

we have checked the stability of our results for different initial values of Zb,k-, h\, and A„ ; a- 
Furthermore, the fixed point solution exists for all considered truncations of the effective 
average potential (see Tab. IIVI) . where we have taken the Taylor series expansion of the 
effective potential up to the tenth order in the field $. The independence of the scaling 
solution both on the initial conditions and higher order operators in the effective potential 
is a manifestation of universality near a continuous phase transition. 

In the symmetric phase both the wavefunction renormalization Zb and the renormalized 
mass m 2 R = Z R l m 2 at the scale k receive large contributions from the massless fermions. 
That is t]b — > 1 for k — > even far from the phase transition which can be clearly seen in Fig. 
[2](b) where the boson anomalous dimension assumes a value close to one in the symmetric 
phase. To compute the critical scaling we introduce the renormalized mass fh\ at a fixed 
scale kc'. 

rh 2 R (k c , Se) = k 2 «(0) - u' kcAcr) (0)) . (80) 
It is given in terms of the first derivatives of the effective average potential in the symmetric 



phase, where the scale k c = r c rrif> is denned via the parameter r c in a standard way [51]. 

The critical exponent v characterizes the divergence of the correlation length at the 
critical point. Here, the correlation length is identified with the inverse renormalized mass 
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4th order 


6th order 


8th order 


10th order 


VB 


0.989 


0.999 


1.003 


1.000 


VF 


0.223 


0.211 


0.207 


0.210 


V 


1.922 


1.936 


1.791 


1.874 


7 


1.942 


1.939 


1.786 


1.875 


1/(2 - n B ) 


1.942 


1.939 


1.786 


1.875 


P 


1.911 


1.935 


1.793 


1.874 


\v(d-2 + r) B ) 


1.911 


1.935 


1.793 


1.874 



TABLE III. N = 2 critical exponents for different orders in the series expansion. 



2ll . |49| as given in (180]) and the critical exponent v is denned as 51| 



- lim 

2 <5e->0 



dhim 2 R {k C} 5e) 



where 



d In 5e 



u(k, Se) 



lim (z>(fc c , Se) + v(k c , Se)u) 

5e->0 



1 d\nm 2 R {k, Se) 



d In Se 



1^1^) 
v ' 1 2 dt 



(81) 

(82) 
(83) 



The value for the critical exponent v is independent of the parameter r c , as long as r c < 1. 
This essentially corresponds to the requirement that the scale k c is sufficiently close to the 
limiting value k — > 0. In our calculations we have taken r c ~ 0.01. 

The critical exponent 7 determines the divergence of the susceptibility which is encoded in 
the non-renormalized mass rh 2 = ZBm 2 R 2l|, |49( . Although it is evaluated in the symmetric 
phase, it is not affected by the fluctuations of the fermions. We have 

<91nm 2 (<5e) 



7 



lim 



(84) 



Se^o d\n5e 

Finally, the critical exponent /3 measures the fluctuations of the renormalized order pa- 
rameter po,fc and is defined in the broken phase: 



1 ^ dlnp 2 



(85) 



2 5e^o dlnSe 

We extract the anomalous dimensions rjB and r/p the same way as the critical couplings. 
Close to the critical parameter e^cr) the renormalization group flow approaches the fixed 
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/I T rl /""\Y*/~i /~\Y* 

^Lll UlutJI 


ULll Ollifcil 


OlLl OI Ciei 


lULIl UIClCl 


e* 


0.4842 


0.5242 


0.5424 


0.5288 




10.7678 


10.3744 


10.1573 


10.3210 






-48.5405 


-73.0962 


-54.6552 


^4,* 






-1956.82 


-485.084 


^5,* 








219713 


hi 


12.8622 


12.9203 


12.9438 


12.9264 



TABLE IV. N = 2 fixed point values for different orders in the series expansion. 

point solution where the system is scale-invariant. That is, the solutions to the flow equation 
stay constant over a wide range of scales where we may extract the corresponding quantities. 
The values of the anomalous dimensions t]b and rjp are defined at the critical point in the 
window where we have a plateau (see Fig. [2j). Our results are summarized in Tab. Ill II where 
we show the values of the critical exponents. They are given for different orders of the 
series expansion for the effective average potential and were obtained using the optimized 
regulator ( |32l) . We have also calculated the critical exponents for the exponential regulator. 
For that calculation, however, we neglect the dependence on the anomalous dimensions in 
the threshold functions. Since the anomalous dimensions are of order one, this gives a very 
rough estimate of the systematic error for our results. We find an agreement of the critical 
exponents on the 10% level. The scaling relations 7 = v (2 — t\b) and (3 = \v {d — 2 + 7]b) for 
the critical exponents also provide a consistency check of our calculations and are also given 
in Tab. IIII1 We see that our results show a reasonable convergence in the series expansion. 

In Tab. II VI we give the values of the critical parameters and couplings, where the asterisk 
denotes the fixed point values, i. e. e* = u'^, A 2i * = etc. These quantities are not universal 
and depend on the particular renormalization group scheme. It is important to comment 
on their behavior in the series expansion of the effective average potential. Taking into 
account only the relevant operators, that is, expanding the potential to fourth order, yields 
a reasonably good result for the scaling exponents. This can also be seen directly in Tab. 
IIVI where the inclusion of higher order irrelevant operators does not significantly alter the 
values for the relevant critical couplings, in contrast to the higher order couplings, that 
vary strongly for different orders of the expansion. The relevant couplings are completely 
stable and show that the important physical information is captured already in the lowest 
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truncation with all relevant operators included. 



VII. CONCLUSIONS 

We have calculated the critical exponents at the quantum critical point for the three- 
dimensional matrix Yukawa type model with U(2) symmetry, which describes N = 2 species 
of Weyl fermions. This can be regarded as a simple low-energy effective model of spinless 
fermions on the honeycomb lattice. The critical exponents of the second order phase transi- 
tion define a new universality class distinct from Gross- Neveu or Neveu- Yukawa type models. 
In particular, this system is special in the sense that it is characterized by large values of 
the anomalous dimensions. Similar results have been obtained in a Dirac cone model for 
graphene where the semimetal-superfluid transition was investigated using functional renor- 



malization group techniques [52|. There, also a second order phase transition was found with 
large values for the anomalous dimensions, both for the anomalous dimensions of composite 
and fermion fields. In the context of compact three-dimensional QED one also observes a 
large value for the anomalous dimension of the gauge field i]a = 1, where the result holds 



exactly due to gauge invariance [53|, |54j. Whether these non-trivial properties can be found 
in suspended graphene is an important open question. To see if these results are indeed 
relevant for graphene requires us to include the long-range Coulomb interactions, which is 
work in progress. 
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Appendix A: Definition of propagators 

In a constant background field Q a b the inverse regularized propagator VZ X is given by 



(z B (l + r B )p 2 + M B ^ 

Z Ftk (l + r F )f -iM F 

Z F (l + r F )i> + iM F y 



(Al) 
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where M B and M F define the scale-dependent mass matrices. 
The background field propagator takes the form: 

( P B {V) ^ 

V k (p)= P { +\p) , 

K P { F )T {p) ) 

where the boson propagator is given by 

P B (p) = (Z B (l + r B )p 2 + M 2 B y 1 , 

and the fermion propagator 

pW(p) = P F (M F ) (Z F (1 + r F )t)TiM F ) 

with 

P F (p) = (Z 2 (l+r F ) 2 p 2 + M 2 y 1 . 



Appendix B: Threshold functions 



For the three-dimensional optimized regulator the shape functions are given by 

r B {y)= Q- lW -y) , 

r F (y) = (^-i) 0{l-y) , 
and the threshold functions can be calculated analytically. They are given by: 

n K ,VB) d V d + 2 J (1 + WY+ 1 



i(F) ( \ 2(5„ i0 + n 



1- VF 



1 i n l VF 



d(l + wi) ni (l + w 2 ) n2 [1 + wi V d + 1 

n 2 ( VB 



m> 2 F \w-,r) F ) 



1 + w 2 V d + 2 
1 



mf\w] T] F ) = — — — + 



1 + w) 4 

1 1 - r) F 1 fl-r}F 



m { J F ^ 2 (w 1 ,w 2 ;r)F,VB) = (l - J^jj 



1 + wY d-2 (l + w) 3 \2d-A A J (1 + w) 2 
Vb \ 1 



d+lj (1 +w 1 ) ni (l + w 2 ) n ^ 
The general form of the threshold functions for generic regulators is given in e. g. 
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Appendix C: Mass eigenvalues 



For the bosonic contribution to the effective average potential Uk we consider only the 
dependence on the quadratic invariant p = |tr$ 2 . The derivatives that appear in the 
evaluation of the mass matrix 



7 



5$ Tab {-p) 



V 



5^ cd {q) 



are given below: 



5<S> Tab (-p) 



$ab(p + V') 



8p(p> 



5 2 p(p'^ 



5$ Tab {-p)5$ cd {q) 



5$ ab (q) 
8ac$bd8(p',-p + q) ■ 



We obtain 



M B m\ cd {p,q) 



5 2 U k 



5p(p') 8p(q') 



x,p,q 



5p(p')5p{q') 5$ Tab {-p) 5§ cd {q) 



+ 



6U k 



8 2 p(p'^ 



x,p> 5p{j>')5<& Tab (-p)5<$> cd (q) ' 



(CI) 



(C2) 
(C3) 



(C4) 



In the non-diagonal constant background configuration $ a b = $oS a b the above expression 
simplifies: 



m 

dp 



T 



S(p,q) . 



(C5) 
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